Excitation of superconducting qubits from hot non-equilibrium quasiparticles 
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Superconducting qubits probe environmental defects such as non-equilibrium quasiparticles, an 
important source of decoherence. We show that "hot" non-equilibrium quasiparticles, with energies 
above the superconducting gap, affect qubits differently from quasiparticles at the gap, implying 
qubits can probe the dynamic quasiparticle energy distribution. For hot quasiparticles, we predict 
a non-neligable increase in the qubit excited state probability P e . By injecting hot quasiparticles 
into a qubit, we experimentally measure an increase of P e in semi-quantitative agreement with the 
model. 
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The superconducting qubit [l|, Q is an excellent candi- 
date for building a quantum computer, as demonstrated 
by recent implementations of the Toffoli gate 0, |3 along 
with simple forms of quantum error correction [5| and 
Shor's algorithm (fjj. However, qubit performance is 
affected by interactions with a range of environmental 
defects, including individual two-level states 0, @| and 
flux noise (ol— Till] . As such, superconducting qubits are 
a sensitive probe of the physics of microscopic defects. 
Non-equilibrium quasiparticles present another impor- 
tant class of defects, which have recently been shown 
to cause qubit decoherence [l2l - fl8j ]. The sensitivity of 
qubits to quasiparticles, and their ability to measure both 
energy emission and absorption rates, enables new mea- 
surements of the non-equilibrium properties of supercon- 
ductors. 

In several recent experiments, the excited state popu- 
lations of superconducting qubits were measured to be in 
excess of thermal equilibrium values, an effect attributed 
to thermal heating induced by non-equilibrium quasi- 
particles [IB-El | . This was supported by the observa- 
tion that the qubit excited-state population was signifi- 
cantly lowered when the level of stray infrared radiation, 
and hence quasiparticle density [llf, was reduced [lij . 
In these experiments, the quasiparticle-induced thermal 
heating necessary to produce the excited state popula- 
tion was thought to result in effective qubit temperatures 
of 70-200 mK [H-0, HI], even though these tempera- 
tures are comparable to the qubit en ergy Eg e ~ 300 mK. 



This violates the typical assumption lfj|, |22|, |24j, |25| that 
[E— A), kgT <C E ge for characteristic quasiparticle ener- 
gies E, superconducting gap A, and dilution refrigerator 
temperature T ~ 20 mK. This means that the specifics of 
the quasiparticle energy distribution cannot be neglected. 

Here, we provide a quantitative theory explaining how 
"hot" quasiparticles can directly excite qubits, provid- 
ing a model for these and other experiments which have 



measured qubits with non-negligible excited state pop- 
ulations [23, 2^-2^] and excitation rates [l^]. We fur- 
ther provide an experimental test of this model by in- 
jecting a non-equilibrium quasiparticle population into 
a superconducting qubit and using the qubit to dynam- 
ically probe this population. We find semi-quantitative 
agreement between our model and the experimental data 
presented here, showing that non-equilibrium quasipar- 
ticles provide a mechanism for the spurious excitation of 
superconducting qubits, distinct from thermal effects. In 
addition, our approach provides a new method to study 
the temporal dynamics of the non-equilibrium quasipar- 
ticle energy distribution and can be used to validate al- 
ternative methods [3fjj . 

A quasiparticle tunneling through the Josephson junc- 
tion barrier in a qubit can cause both excitation and dis- 
sipation in the qubit, as illustrated in Fig.[TJ Consider a 
qubit initially in its excited state: A "cold" quasiparticle 
near the gap energy can absorb the qubit transition en- 
ergy E ge between the qubit's excited |e) and its ground 
\g) states, causing the qubit to switch to its ground state 
(blue arrows in Fig.[Tfa) and (b)). Any quasiparticle in 
the junction area can absorb this energy, so the qubit 
|e) — > I g) decay rate T± due to this channel is proportional 
to the quasiparticle density n qp . For a qubit initially in 
its ground state, a "hot" quasiparticle sufficiently above 
the gap energy can excite the qubit, but only if the quasi- 
particle has energy greater than A + E ge (red arrows in 
Fig.QJa) and (c)). The qubit \g) — > |e) excitation rate 
thus depends on the energy distribution of the quasipar- 
ticle population. 

If the quasiparticle population were well-described by 
a temperature T ~ 20 mK -C E ge /kB, then a negligi- 
ble qubit excitation rate T-f- would be expected, as in 
Fig.QJb). There are however a number of processes that 
can produce quasiparticles with energies well above kpT, 
which then relax via quasiparticlc-phonon scattering [31| , 
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FIG. 1: (Color online) Qubit decay F^ (blue) and excitation 
Tf (red) mediated by tunneling quasiparticles. (a) Portion of 
the superconducting phase qubit potential energy diagram, 
showing the F^ and F^ transitions between the ground state 
\g) and the excited state |e), along with the qubit energy E ge . 
(b) Cold non-equilibrium quasiparticles, which have energies 
near the superconducting gap A, can only absorb E ge , result- 
ing in qubit F^ decay. The density of states (p(E), horizontal) 
on both sides of a Josephson junction (superconductor S - in- 
sulator I - superconductor S) is shown versus quasiparticle en- 
ergy E (vertical). The quasiparticle energy distribution f(E) 
is shown by the shaded triangles, (c) Hot non-equilibrium 
quasiparticles with energy above A + E ge (portion of f(E) 
with E > A + E ge ) not only can cause qubit transitions 
but can also relax by causing qubit \g) |e) transitions. 



ing as the square root of the injection rate; we verify be- 
low (Fig. [3]) that the qubit excited state probability P e 
is independent of the injection energy. Steady state is 
achieved by balancing phonon scattering and quasiparti- 
cle injection and recombination [32|. Although the result- 
ing steady-state occupation probability f(E) cannot be 
described by an effective temperature, its dependence on 
quasiparticle energy for A < E < 1.4 A is similar [33| to a 
thermal distribution with temperature of approximately 
70 mK. 

With f(E) determined in this manner, we calculate 
the qubit excitation rate and decay rate For 
a tunnel junction with resistance Rt and capacitance 
C, and for normalized quasiparticle density of states 
p(E) = E/y/E 2 — A 2 , the qubit decay rate induced by 
all quasiparticles is [FJ, [14 1 



r± = 



RtC 



dE EEf + A 5 



E, 



y>- 



EE f 



p(E)p(E f )f(E), 



(1) 

where the final quasiparticle energy Ef = E + E ge is 
higher than E due to the absorption of the qubit energy 




but for which the non-equilibrium quasiparticle occupa- 
tion probability f(E) still has a significant population 
of hot quasiparticles [F3| , with energies well above fc^T. 
The distribution f(E) is typically not thermal. 

In order to model the steady-state quasiparticle dis- 
tribution, we assume quasiparticles are injected in the 
junction at a constant rate at an energy Ei n j well above 
A + E ge , with the resulting quasiparticle density n qp seal- E, 



qpl cp 

FIG. 2: (Color online) Calculated effects of hot quasiparticles. 
Quasiparticle densities incorporate scattering and recombina- 
tion, with quasiparticles injected at energy Einj = 1.8 A (solid 
lines) or A + 0.94_E 9e (dashed line), assuming A/fcs = 2K as 
for aluminum. Plots are for three different qubit frequen- 
cies E ge /h: 13.3 GHz (black), 6.7 GHz (green), and 3.3 GHz 
(red); lines are guides to the eye. (a) Qubit decay rate Tj. 
normalized by the cold rate FJ from Eq. ([2]), showing that TJ 
is a good approximation to T^. (b) Qubit \e) state probabil- 
ity P e — r-|-/(rj, + r-f) vs. normalized quasiparticle density 



o/n C p. The je) state occupation increases linearly with the 
quasiparticle density for low densities. Injecting low energy 
quasiparticles (dashed line) results in a greatly reduced P e . 



ge . Here, cj) is the junction phase, which is typically 
4> ~ for the transmon and <f> ~ ir/2 for the phase qubit. 
For cold non-equilibrium quasiparticles, corresponding to 
f(E) having population only at the gap energy A, this 
integral gives 



1 + cos 4> i A 
V2 R T C^Eg~e 



3/2; 



(2) 



Here, n cp = D(Ep)A is the Cooper pair density, 
n qp = 2D(E F ) p(E)f(E)dE is the quasiparticle den- 
sity, incorporating both hot and cold quasiparticles, and 
D{Ep)/2 is the single spin density of states. This is 
the standard result for quasiparticle dissipation 11511 . and 
is equivalent 14{ to the Mattis-Bardeen theory [34j for 
= 0. 

In Fig.rjja) we plot T^/T9, the ratio of the numerically- 
integrated qubit decay rate T± assuming both hot and 
cold non-equilibrium quasiparticles (Eq.[l} to the rate Pr 
for cold quasiparticles at the gap (Eq.[2]), with f{E) cal- 
culated as explained above. We see that Tj, s» Pf for a 
range of parameters, so the quasiparticle-induced decay 
rate is determined primarily by the total quasiparticle 
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FIG. 3: Calculated qubit excitation probability P e vs. quasi- 
particle injection energy Ei„j. The injection energy is ex- 
pressed as the fraction of the qubit energy E ge = h(6.7 GHz) 
above the superconducting gap A. For Emj > A + 1.7E ge , 
P e is essentially constant. Here, the injected quasiparticle 
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density n qp and depends only weakly on the quasiparti- 
cle occupation distribution. 

However, the quasiparticle distribution is key to char- 
acterizing the quasiparticle-induced steady-state excited 
state population, P e — T-^/iT^ + T^). The rate can 
be evaluated using Eq. (pQ) by substituting Ej = E— E ge , 
changing the lower limit of integration to A + E ge , and 
dropping the factor (1 — f(Ef)) as it has negligible effect. 
Using the same f(E) as before, we calculate the proba- 
bilities plotted in Fig.[2jb). Notice that a non- negligible 
probability of a few percent can be obtained for quite 
modest quasiparticle densities. The probability P e de- 
creases for smaller quasiparticle densities and for larger 
qubit energies, as expected. For small occupation prob- 
abilities, this result can be approximated by the fit func- 
tion P e ~ 2.17{n qp /n cp )(E ge /A)- 3 - e5 . 

To determine the sensitivity of this result to the quasi- 
particle injection energy, we also calculated P e as a func- 
tion of the injection energy Ei n j. As shown in Fig.[3j 
P e is independent of the injection energy for Ei n j > 
A + 1.7E ge ; hence, for sufficiently large injection ener- 
gies, the actual injection value is unimportant. In addi- 
tion, we see that P e is significantly suppressed for quasi- 
particle energies below A + E ge , demonstrating that cold 
quasiparticles do not excite the qubit. The maximum 
at Ei n j = A + E ge is caused by the peaked final state 
density of states p(Ef) at Ef = A in the expression for 

r t . 

To experimentally test these concepts, we performed 
an experiment in which we deliberately injected quasi- 
particles into a phase qubit and measured the excited 
state probability P e and the increase in the qubit decay 
rate ST^, which is proportional to n qp /n cp (Eq.[5]). As 
shown in Fig. [4] we generated quasiparticles by apply- 
ing a voltage pulse above the gap voltage A/e to the 
qubit 's readout superconducting quantum interference 
device (SQUID), thus changing the quasiparticle density 



and f(E). This is similar to previous work DJJ, except 



here the qubit measurements include the enhancement 



FIG. 4: (Color online) Experimental apparatus and proto- 
col, (a) Photo of phase qubit, with voltage V 3q applied to the 
readout SQUID and flux bias <&bias- Two paths for quasipar- 
ticles to diffuse to the qubit junction are indicated by dotted 
and solid arrows, (b) Pulse sequence. Quasiparticles are gen- 
erated with a voltage pulse on the SQUID line of amplitude 
Vinj > 2(A + E ge )/e for a time tinj, which are varied to ad- 
just the quasiparticle density. After a delay time tdif, 40/zs 
of which is at the operating bias, allowing quasiparticles to 
diffuse to the qubit, the qubit state is projected with a flux 
bias pulse and read out using a voltage pulse on the SQUID. 



5P e as a function of the energy relaxation rate increase 
ST i. As shown in Fig.@|b), we first applied a voltage 
pulse to the SQUID to generate above-gap quasiparticles; 
the duration ti n j of this pulse was varied to adjust the 
quasiparticle density. We reset the qubit into the \g) state 
using &bi as and then waited a variable time tdif follow- 
ing the pulse, giving the quasiparticles time to diffuse to 
the qubit and allowing study of the temporal dynamics. 
After the qubit control pulses, we measured the qubit P e , 
reading out the qubit by increasing V sq to approximately 
0.7 'A/e to switch the SQUID into the normal state while 
minimizing quasiparticle generation. 

In Fig. [SI we plot the observed changes SP e versus ST^ 
for two different phase qubits d, |35| as we varied tdif 
and ti n j over the ranges indicated in Table I. We find 
that SP e monotonically increases with 6T±, thus scaling 
with quasiparticle density, as predicted by the theory. In 
fact, by deliberately injecting hot quasiparticles, we were 
able to increase P e by more than 10%; this directly shows 
that hot quasiparticles can significantly excite the qubit. 

The quasiparticle density was changed in three ways: 
Varying the diffusion time tdif (triangles), the injection 
time t^j, and the injection voltage Vinj (open vs. closed). 
There is essentially no difference between changing the 
injection voltage and the injection time; this makes sense, 
as P e is relatively insensitive to the injection energy 
(Fig. [3]). However, we observed a significant difference 
between varying the diffusion time and varying the injec- 
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FIG. 5: (Color online) Experimental impacts of injected 
quasiparticles. Increases in qubit decay rate <5P^ (proportional 
to quasiparticle density n qp /n cp ) and the raw qubit excited 
state probability 5Pe are measured with respect to the values 
without quasiparticle injection. The quasiparticle density was 
parametrically varied by changing either the length tmj of the 
injection pulse (circles) or the diffusion time tdif (triangles); 
arrows indicate ti n j + tdif ~ 400 pis. Data are for two slightly 
different designs, as described by (a) [(| and (b) [3f|; symbols 
and experimental parameters are defined in Table Theoret- 
ical predictions (solid lines) were calculated for quasiparticles 
injected at I.9A using measured E ge , giving E ge /A ~ 0.15. 



tion time. This is not surprising, as P e is sensitive to the 
quasiparticle energy distribution, which changes with the 
diffusion time. As shown in Fig.0^, quasiparticles gen- 
erated at the measurement SQUID must diffuse through 
approximately 700 /im of ground plane metal in order to 
reach the qubit junction. This gives time for quasiparti- 
cle scattering and relaxation with respect to the injected 
distribution. In fact, those triangles in Fig. [5] where 5P e 
begins to level off correspond to Unj + t<u / ~ 400 fjs (de- 
noted by arrows). This is of the same order as typical 
quasiparticle recombination times 3^, 37 1. where f(E) is 
expected to become dominated by thermal equilibrium 
excitation independent of quasiparticle density. However, 
for times less than this (see Table H] for times), the data 
where tdif was varied is more similar to the data where 
Unj is varied, indicating a non-thermal f(E) dependence 
of P e which varies with time. 

Predictions from the model are plotted in Fig.[5]as solid 
lines. Although the power-law dependence (slopes) be- 
tween model and measurement are in reasonable agree- 
ment, the model 5P e is low by about a factor of three 
for qubit (a) and high by about a factor of four for qubit 
(b). This can not be attributed to the small differences in 
the qubit geometries, because a third device (not shown), 
with the same design as (a), gave data about a factor of 
two lower than the model. These differences could be 
due to subtleties in the model not considered here. For 



TABLE I: Parameters for data in Fig. [5] Here we list the pa- 
rameter, tdi f or ti n j , that is varied for the five data sets, along 
with its values from the left to the right sides of the figure. 
The value of the fixed parameter and the injection voltage 
Vinj are also provided. Experimental parameters measured 
independently without quasiparticle injection include P e , T^, 
qubit frequency E ge /h, and qubit critical current 7 C . 
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instance, there can be sample-to-sample variation such 
as a difference in film thickness or film quality, yielding 
variations in the gap energy, in the quasiparticle diffu- 
sion path, and in junction parameters. In addition, we 
note that the calculation of P e assumes non-equilibrium 
quasiparticle relaxation through electron-phonon scat- 
tering and recombination. However, other effects such 
as quasiparticle diffusion and trapping of quasiparticles 
from non-uniform gaps may significantly affect the dis- 
tribution f(E), altering the prediction for the excita- 
tion rate. In addition, even though quasiparticles are 
expected to be generated with energies of eVi n j/2, in re- 
ality there will be a distribution of quasiparticle energies 
centered around this value for a given Vinj, so even for 
eVinj > 2(A + Ege), some of the generated quasiparti- 
cles may instead be cold. We conclude that the simple 
theoretical model used here for the quasiparticle energy 
distribution f(E) is only semiquantitative, predicting 
in the non-thermal regime to about a factor of 4. 

In conclusion, we have shown that "hot" quasiparti- 
cles with energies greater than A + E ge can cause sig- 
nificant ground-to-excited state transitions in supercon- 
ducting qubits. This is in contrast to "cold" quasipar- 
ticles solely at the gap, which only cause superconduct- 
ing qubits to relax. This means quasiparticles cannot 
be adequately described by a single parameter such as 
n q p/n cp or temperature. As illustrated by varying both 
U n j and tdi f , the particular quasiparticle distribution af- 
fects the observed qubit excitation probability. This the- 
ory semiquantitatively matches the observed behavior of 
the qubit excitation probability versus quasiparticle den- 
sity. 
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